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Abstract 

The validity of small damping approximation (SDA) for the quasi-classical description of the 
averaged properties of nuclei at high temperatures is studied within the framework of collisional 
kinetic theory. The isoscalar collective quadrupole vibrations in hot nuclei are considered. We show 
that the extension of the SDA, by accounting for the damping of the distribution function 5f in 
the collision integral reduces the rate of variation with temperature of the Fermi surface distortion 
effects. The damping of the 5f in the collision integral increases significantly the collisional width 
of the giant quadrupole resonance (GQR) for small enough values of the relaxation time. The 
temperature dependence of the eigenenergy of the GQR becomes much more weaker than in the 
corresponding SDA case. 

PACS numbers: 21.60.Ev, 24.30.Cz 
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I. INTRODUCTION 



In this work, we study the damping of nuclear multipole vibrations within the semi- 
classical kinetic theory. Semiclassical approaches are quite instructive for the description 
of averaged properties and macroscopic collective motion in multiparticle systems, see Ref. 

and references therein. In many cases, they allow us to obtain analytical results and 
represent them in a transparent way. However a significant simplification is achieved in 
these approaches due to the so-called small damping approximation (SDA), which assumes 
that the relaxation rate is small enough compared to the energy of the collective excitation. 
This approximation is reasonable in a cold Fermi system where the excitation energy is high 
enough (due to the distortion of the Fermi surface) and the relaxation rate is quite small 
(zero sound regime) ODD . With the heating of the Fermi system the energy of the collec- 
tive vibrations decreases and the two-body relaxation rate increases simultaneously (zero- 
to first- sound transition) Q|. Thus with increasing temperature the ratio of the relaxation 
rate to the excitation energy increases rapidly and the use of the SDA can be invalid. In 
this case the analytical results for eigenenergy, width, strength function, etc., obtained with 
the SDA, have to be verified independently. 

In this paper we study the validity of the SDA within the Fermi liquid theory considering 
the excitation of the giant quadrupole resonance (GQR) in a hot nucleus. In Sect. II 
we derive the equations of motion for the collective quadrupole vibrations in a finite heated 
Fermi liquid drop for both the SDA and the general case. In Sect. Ill we provide a numerical 
study of the temperature dependence of the excitation energy and the relaxation rate for 
the above mentioned cases. Our conclusions are given in Sect. IV. 



II. EQUATIONS OF MOTION 

To derive the equations of motion for the nuclear shape variables, we start from the 
following kinetic equation for the small deviation Sf of the distribution function / = /(r, p; t) 
from the one in equilibrium, / eq , taking into account the collision integral SSt linearized with 
respect to equilibrium jj, [J , 

pf + L6f = 6St[5f]. (1) 
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On the left hand side of Eq. (|T|). the operator L represents the drift terms including the 
selfconsistent mean field V, 

L5f = ^ • V r 5f - V r V cq ■ V p 5f - V r 6V ■ V p / eq • 

m 

We will follow the nuclear fluid dynamic approach, see Ref. [J, and take into account the 
dynamic Fermi surface distortion up to multipolarity 1 = 2, 

1=2 

6f = ~(%) ^f^m^ip). (2) 

Here e is the quasiparticle energy Q]. A generalization of this approach to the case of an 
arbitrary multipolarity / of the Fermi surface distortion can be done in a straightforward 
way, see Ref. p. Using Eqs. (JTJ) and (J2J) we will derive a closed set of equations for the 
following p-moments of the distribution function, namely, local particle density p, velocity 
field u v and pressure tensor P u ^, in the form (for details, see Refs. jilB, Isj]) of 

<><)p + 4t(p^) = °> ( 3 ) 



dt dr 



V 



+ + (4) 



dt q <9r Q \5p 2 J eq dr v 



9P Lp . D / du a dup 2 du u 



Here £ is the internal energy density, which is the sum of the kinetic energy density of 
the Fermi motion and the potential energy density associated with the nucleon-nucleon 
interaction. The equilibrium pressure of a Fermi gas, P eq , is given by 



3m J (2tt^) 3 

P' a p is the deviation of the pressure tensor from its isotropic part due to the Fermi surface 
distortion 

Kp = ~~ / (2^) 3 ^° ~ ~ mu ^ 5Z 5 /2m(r,t)r 2m (]3), (7) 

SSta/3 is the second moment of the collision integral 

6Sta(3 = m" / J0f P*Pd 6St i 6 f] ( 8 ) 



and g = 4 is the spin-isospin degeneracy factor. In Eqs. ©-© and in the following 
expressions repeated Greek indices are to be understood as summed over. 
Using the Fourier transformation for the pressure 

PUt) = I %^ P ^ (9) 

and similarly for the other time dependent variables, we find the solution to Eq. (j5J) as 



P'al3,u> — Y^J-— :-PeqA a/ 3 itlJ , (10) 



where we used the symbol 
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for this combination of gradients of the Fourier transform Xa,w of the displacement field. 
The time derivative of the displacement field x{ r >t) is defined as the velocity field, hence 

U a ,uj = -iux<x,u- (12) 

To obtain Eq. (fTU|) we have also used the fact that the tensor (JSt^, Eq. (jHJ), can be reduced 
to 

SStapu = P a B.wi (13) 

T 

due to our restriction to quadrupole deformation of the Fermi surface. This is because the 
I = and 1 components of the expansion (J2J) do not contribute to the collision integral, 
reflecting the conservation of particle number and momentum in a collision, respectively. 
The relaxation time r in Eq. (|13|) is derived by [t| 

1 [dp p 2 Y 20 (p)6St[5f] 

fdppiY 20 (p)5f ■ 1 ) 

The collision integral 5St[5/] is related to real transitions in the interparticle collisions 
and the expression for this collision integral should contain ^-function associated with the 
conservation of energy. Note that the form of Eq. (TEH) is also correct for a non-Markovian 

n 

collision term, with the collision time r being dependent on the frequency u, see Ref. j8|. 
This ^-dependence is associated with the memory effects in the interparticle collisions and 
the collision integral takes then the following form p] 



SSt[5f] 



2(27rfi)« W{iPj}) 5iAp %W> 



5fj [5(Ae + hu) + 5(Ae - hcu)] . (15) 

eq 



Here, ft = /(r,p;;t), 



Ae = ei + e 2 - e 3 - e 4 , Ap = pi + p 2 - p 3 - Pa, 

where ej = p 2 /2m + V eq (r) is the classical single-particle energy, Q is the Pauli blocking 
factor 

Q = (l - - f 2 )hft - A/ 2 (i - / 3 )(i - /«), 

^({Pj}) = {dtJin/ dQ)4g(2irh) 3 /m 2 is the probability of the scattering of nucleons near the 
Fermi surface and da in /dVL is the corresponding in-medium differential cross section. 

The collision integral of Eq. (I15 | ) provides the relaxation time tsda of small damping 



approximation. It is given by 



3, LLC 



«5 rifi^i 



where T is the temperature and a and ( are parameters. The phenomenological parameter 
a is determined by the in-medium spin-isospin averaged nucleon-nucleon cross section o~ in . 
The magnitude of parameter ( = 1 was earlier obtained by Landau [Uj in the calculation 
of the absorption coefficient in a Fermi liquid. The parameter ( can be also evaluated in a 
direct way within the framework of a kinetic theory if the exact expression for the collision 
integral is known. In a general case, the tu-dependence (memory effect) in JSt[5/] is caused 
by three sources |7j: (i) variation of the distribution function 5f, (ii) variation of the mean 
field SV, and (in) screening effect for the two-body scattering in hot Fermi system due to 
high frequency collective vibrations. For the sake of simplicity we will only consider the 
contribution to 5St[5/] from the variation of the distribution function 5f. The constant £ 
is then given by ( = 3 0. 

The frequency u in Eq. (fTBj) is assumed to be real. This reflects the fact that we 
neglect the damping of the distribution function 5f in Eq. ()15|) due to the small damping 
approximation. In the case of the complex eigenfrequency to = ujr — iui with non-vanishing 
imaginary part uj > the delta function in Eq. (fTKj) should be spread. Thus, to extend the 
SDA by taking into account the damping of the distribution function 5f , it is reasonable 
to replace the ^-functions in Eq. (|15|) by Lorentzians of the form 

6(Ae±hu)=>- T - A kU ' - (17) 
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Using Eqs. (JH|), (fH>j) and (JT7J), one then obtains, see also Ref. [12 

4:7T 2 ah 



T 



RW + RU 



de(e-V cq f 2 (df/de 



Icq 



(18) 



where 



(±) _ _3hujj_ f°° de 1 de 2 de 3 de i ( £l - V cq f' 2 (h5_Q_dh \5_Q_dh 5Q_dj\ 
' ' 2tt L (Ae ± huj R ) 2 + (hujj) 2 \2 5h de x 2 5f 2 de 2 5f 3 de ; 



(19) 



cq 



Here and below we assume a constant spherical potential of radius Rq for V eq and use the 
Fermi function for the equilibrium distribution function / eq 

' p 2 /2m — A x 1 



cq 



1 + exp 



T 



(20) 



A ~ e F 



12 W 



where p 2 /2m = e — V eq , e F = p 2 F /2m and pf is the Fermi momentum 

We will reduce the local equations ©-© to macroscopic equations of motion for the 
nuclear shape variables assuming the following separable form for the displacement field 

X(r,*) = /3(*)v(r), 

where (3{t) = @ u exp(-iujt) is the collective variable. Considering quadrupole vibrations, we 
will assume an irrotational motion with 



v(r) = V(r 2 F 20 (f))/2. 



(21) 



Substituting Eq. (|2*T|) into Eq. (jlj) we obtain the equation of motion for the collective 
variable (3(t), 

mp^vj+lp^l 6 ^] S Pw + d -^](3 = 0, (22) 



dr a \5p 



cq 



dr h 



where 



Multiplying Eq. (|2*2*|) by v a , summing over a, and integrating over r-space, we obtain the 
dispersion equation for the eigenfrequency u, 



-Buj 2 + C LD + C'{u) - iuA{u) = 0. 



(24) 
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Here, B is the mass coefficient with respect to the collective variable /3(t): 

B = m J dr p cq (r)|v(r)| 2 . (25) 

The dissipative term A(u) and the stiffness coefficient C'(u) are caused by the dynamic 
distortion of the Fermi surface and they are given by 

A{lo) = [ dr — Im n uattl3 V u v a , (26) 

C'(u) = J dr ( ^- Im n ua ^ - Re Tr vait A V v v a . (27) 

The stiffness coefficient Cld in Eq. (J23|) can be identified with that of the traditional 
liquid drop model [3| 

where bs = 17.2 MeV and be = 0.7 MeV are the surface tension coefficient and the Coulomb 
energy coefficient, respectively. To evaluate the kinetic coefficients A(u), B and C'(u) we 
will assume a sharp surface of the nucleus in equilibrium p eq = p 9(R — r), where p is 
the bulk density and Rq = r A 1//3 fm. Taking into account the expression (|21|) for the 
displacement field, we obtain from Eq. (J25|) the mass coefficient as 

B = —AmRl. (29) 
8n 

Using Eqs. ©, JED, (HU), (ED) and fl23J), we obtain from Eqs. flUJ) and (J27D the kinetic 
coefficients A(uj) and C'(oS) in the following form 

= (i-^)Wf J ' (30) 



Here, 



and 



J = J drP cq A aP Vpv a = 5R 3 P (32) 



2 /«_2\ 2/3 



lr /67T 



P = 3- — P 

am \ g 



5/3 
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III. NUMERICAL CALCULATIONS 



For the numerical calculations in this work we adopt the value of tq = 1.12 fm. In Fig. 
1, we have plotted the results of the calculations of the excitation energy E = hu^ and 
the collisional width T = 2hujj obtained from Eq. (J24)) for the nucleus with A = 224 and 



Z = A[l -2A/5(A + 200)]/2 = 88, which corresponds to the valley of beta-stability 13]. We 



have here adopted the value of a = 9.2 MeV from [7], which corresponds to <j- m « a frce /2, 
where <7f ree ~ 40 mb is the cross section for the nucleon-nucleon scattering in free space. The 
eigenenergy E in a cold nucleus is strongly shifted up with respect to the one in the nuclear 
liquid drop model because of the contribution from the term C'(u) to the nuclear stiffness 
coefficient C = Cld + C'{u) in Eq. (|24|). The stiffness coefficient C'(u) is caused by the 
dynamic Fermi surface distortion effect and it decreases strongly with increasing temperature 
T. Due to this fact, the stiffness coefficient C and the eigenenergy E decrease monotonously 
with temperature and approach the liquid drop model limit for large temperatures, see also 
Ref. 0. As seen from Fig. 1, the SDA leads to a significantly faster decrease of E with 
increasing temperature T, i.e., the SDA overestimates the decrease rate of the Fermi surface 
distortion effect with temperature. The damping of 5f in the collision integral does not 
play an essential role in the derivation of the collisional width T of the GQR. We can see 
from Fig. 1 that the damping of Sf increases the width T with respect to the SDA, but the 
additional contribution to the width does not exceed 20%. 

The difference between the SDA and its extension of Eq. fjlSjl becomes stronger for 
smaller values of the parameter a. Let us introduce the ratio 

<? = tsda/V, (33) 

which allows one to estimate the value of the correction for the two-body relaxation time 
due to the damping of the distribution function Sf in the collision integral. The result of 
the calculations of q as a function of the temperature T, with ujr and ujj taken from Fig. 1, 
is shown in Fig. 2 for two values of parameter a = 9.2 MeV and a = 4.6 MeV. (Note that 
the relaxation time r is a complicated function of a because the quantities ujr and uj in 
Eq. (fTU|) depend on a due to Eqs. (|24|) . (|30j) and (|31|h) In Fig. 3 we show the temperature 
dependence of the ratio 

Ar/r = (r - r SD A)/r, 
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for two values of the parameter a. Here T = hojj is the exact collisional width obtained 
using the relaxation time r of Eq. (JTSJl and Tsda is the one obtained within the SDA using 
tsda from Eq. The non-monotonic behavior of the ratio AT/T in the case of a = 4.6 

MeV (dashed line) is because the width Tsda has a bell-shaped form as a function of T, for 
moderate temperatures |?|. 



IV. CONCLUSIONS 



In this paper we have studied the validity of the small damping approximation for the 
quasi-classical description of the averaged properties of nuclei at high temperatures in the 
framework of the Fermi liquid theory. The basic equation of the collisional kinetic theory 
has been reduced to macroscopic equations of motion for the shape variable of the nucleus. 
We have extended the small damping approximation taking into account the damping of the 
distribution function 5f in the collision integral <5St[5/] and replacing the energy conservation 
^-functions in <5St[5/] by Lorentzians. Such kind of an extension of the SDA leads to an 
increase in of the relaxation rate. The corresponding decrease of the collisional relaxation 
time r does not exceed 20% of the value tsda obtained within the SDA and it decreases 
with increasing temperature T. 

We have shown that the extension of the SDA due to the damping of 5f in the collision 
integral reduces the rate of variation with temperature of the Fermi surface distortion effects. 
The smoothing of the ^-function in the collision integral of Eq. (fTo'j) caused by the damping 
of the distribution function 5f provides much more weaker dependence of eigenenergy E on 
the temperature than the one in the SDA. 

The difference between the SDA and its extension of Eq. (fTHj) becomes stronger for 
smaller values of the relaxation parameter a. The dashed line in Fig. 3 shows that the 
correction to the SDA due to the damping of the Sf in the collision integral can increase 
the collisional width T of the giant quadrupole resonance for a small enough value of the 
parameter a. 
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Figure captions 



Fig. 1. Temperature dependence of the collisional width T = 2fiuJi and excitation energy 
E = TujOr of the GQR for the nucleus with mass number A = 224 obtained from Eq. (|2"l*|) 
with a = 9.2 MeV. The solid lines are obtained using r from Eq. (|18|) and the dashed lines 
are for the SDA with Eq. (JTHJ). 

Fig. 2. The dependence of the ratio q for the GQR for the nucleus with mass number 
A = 224 on the temperature T obtained from the solution of Eq. (j24j) . The solid and 
dashed lines correspond to the calculation with the parameter a = 9.2 MeV and a = 4.6 
MeV, respectively. 

Fig. 3. Temperature dependence of the relative correction AT /T to the width of the 
GQR due to the damping effect in the collision integral for the nucleus with mass number 
A = 224. The solid and dashed lines correspond to the calculation with the parameter 
a = 9.2 MeV and a = 4.6 MeV, respectively. 
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